Electron-Mediated Nuclear-Spin Interactions Between Distant NV Centers 
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We propose a scheme enabling controlled quantum coherent interactions between separated nitrogen- vacancy 
centers in diamond in the presence of strong magnetic fluctuations. The proposed scheme couples nuclear qubits 
employing the magnetic dipole-dipole interaction between the electron spins and, crucially, benefits from the 
suppression of the effect of environmental magnetic field fluctuations thanks to a strong microwave driving. This 
scheme provides a basic building block for a full-scale quantum information processor or quantum simulator 
based on solid-state technology. 
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The spins of single dopants in solids are key elements in 
the development of solid-state quantum-information technolo- 
gies [1,2]. In particular, nitrogen-vacancy (NV) colour cen- 
ters in diamond are promising quantum processors: single de- 
fects can be detected using confocal microscopy [3, 4], their 
spin state can be initialized, manipulated, and readout opti- 
cally [5-8], and their quantum coherence survives at room 
temperatures [9] . One of the remaining challenges is to con- 
trol the spin-spin interactions to perform quantum-logic op- 
erations, and major steps along this direction have already 
been accomplished. The hyperfine coupling between the NV 
electron spin and the nuclear spins of neighboring impurities 
( 13 C, 15 N) offers a unique opportunity to build small quantum 
registers [8, 10, 12-14]. These devices can be scaled up by 
means of ion implantation techniques, yielding periodic ar- 
rays of NV centers [15]. However, the controlled couplings 
now require longer-range interactions, as provided by optical 
channels [16], or magnetic dipole-dipole couplings between 
the electron spins [17]. 

Although the feasibility of the magnetic-coupling approach 
has been demonstrated recently [17], fabricated NV arrays 
often suffer from shorter electron coherence times that af- 
fect the fidelity of the quantum gates. From this perspec- 
tive, 14 N or 15 N nuclear spins would be better-suited qubits 
due to their longer coherence times, together with the avail- 
ability of single-shot readout [10]. Unfortunately, the direct 
nuclear dipole-dipole interaction is negligible, which necessi- 
tates the search for alternative schemes to couple the nuclear 
spins. This letter presents a theoretical proposal for imple- 
menting robust quantum gates between two distant nuclear- 
spin qubits mediated by the long-range dipolar interaction be- 
tween electron spins. The main idea is to exploit the long 
nuclear coherence times for storage, and to use the electronic 
degrees of freedom as a quantum bus that mediates the nu- 
clear spin interaction. Such a general scheme can be applied 
to different setups, and has also been proposed for quantum- 
Hall systems [11]. Active control of the spins via microwave 
fields allows reaching high fidelities, even in the presence of 
the magnetic noise associated to the complex mesoscopic en- 
vironment of solid-state systems. In fact, the nuclear driving 
acts as a continuous decoupling mechanism [18] that mini- 
mizes the effects of the noise, and provides a new tool in ad- 
dition to pulsed techniques [19]. 

The model- We consider two NV defects j — 1,2, whose 
unpaired electrons form a spin-triplet ground state Sj = 1, and 



focus on 14 N with a nuclear spin Ij — 1. The Hamiltonian that 



/(e-n) 



describes each NV center is Hj = Hf +HJ l) +H 
Hf =D j ((S)) 2 -^)+ ge iiBB-S j , 
-Pj({I)) 2 -\l^-g n ^B-I h 



H 



(n) 



7 (e-n) 



l -I 2 



(1) 



Hr=AtfI^±Af(Sp 7 +SjIf) 



where S/,// are the electronic and nuclear spin-1 operators, 
and Sf = S x - ± iS y -, 1^ = E ± U y the usual ladder operators. 
Here, Dj(Pj) stands for the zero-field splitting of the electronic 
(nuclear) ground state, B is an external magnetic field, /Ib (Mn) 
is the Bohr (nuclear) magneton, and g e (gn) is the electron (nu- 
clear) g-factor. The electron-nuclei interaction is quantified 
by the hyperfine longitudinal (transverse) coupling A li - (Aj~). 
The present discussion is focused on a single pair of closely- 
spaced NV centers, and we use the realistic parameters of the 
experiment in [17]. We emphasize, however, that this scheme 
can be extended to arrays of implanted NV centers, provided 
that their distance is small enough. Let us also remark the hi- 
erarchy of couplings, Dj > Pj > A U j,AJ, and g e MB > gnMN 
(see Table I, where h = 1). Finally, we introduce the secular 
dipole-dipole interaction between the electron spins 
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where Jn — §Ii^b(^ — 3 cos 2 Q\2)/2cr\ 2 in gaussian units, r\2 
is the distance between the NV centers, cos B\2 = e z • r^/ru, 
and c is the speed of light. For the distances reached in the ex- 
periment, 7*12 ~ 10nm, the dipolar coupling J\2 ~ 70kHz is the 
smaller energy scale in the problem. As mentioned above, the 
magnetic dipole-dipole interaction between the nuclear spins 
is completely negligible since (gnAWgeMB) 2 ~ 10 -8 , and an 
indirect mechanism for the nuclear coupling is thus required. 
Effective static interactions.- In Fig. 1(a), we represent 
schematically the process leading to nuclear spin-spin inter- 
actions. The hyperfine interaction couples the nuclear to the 
electronic spins of each NV center, which are in turn cou- 
pled through the magnetic dipole-dipole interaction. There- 
fore, one may use the electrons as a bus to mediate the 
nuclear coupling. A naive estimate of this coupling fol- 
lows from Fig. 1(b), where we represent the energy spec- 
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trum of Ho = Zj(H) + H j) + #12 • Due t0 the energy- 
scale hierarchy in Table I, the levels are clustered in mani- 
folds determined by the electronic spins |rai,W2) e . The dy- 
namics within the ground-state manifold, |0,0) e , corresponds 
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Figure 1 . Effective nuclear spin-spin interaction, (a) Schematic 
diagram of the electron-mediated interaction between the nuclear 
spins, which exploits the magnetic dipolar interaction and the local 
hyperfine coupling, (b) Diagram of the energy levels of the two NV 
centers. Since Dj is the largest energy scale, the energies are clus- 
tered in manifolds determined by the electronic spin. The transverse 
part of the hyperfine coupling Aj- induces transitions between dif- 
ferent manifolds, and mediates an effective XX interaction between 
the nuclear spins, (c) Schematic diagram of the Zeeman splitting 
for the electronic and nuclear energy levels. By carefully selecting 
the microwave frequencies, we drive a particular electronic and nu- 
clear transition, (d) Energy levels of the driven Hamiltonian. For 
very strong driving £2 e , the electronic spin in the lowest manifold is 
|— ) oc |0) — | — 1). The hyperfine coupling A 1 '- induces virtual transi- 
tions to the excited manifold, split by the dipolar interaction and the 
inhomogeneous broadening, and leads to an effective ZZ interaction. 



to nuclear spin flips |Mi,M 2 ) n -> \M' v M' 2 ) n , with Mj,Mj = 
0, ±1, and follows from second-order processes where the hy- 
perfine coupling virtually populates states from the excited 
manifold. Therefore, a crude estimate of the dynamics is 
#eff ~ Jefttfl^ + H.c, where J G ^ oc (A^A^/D. A more care- 
ful Schrieffer- Wolff-type calculation takes into account the 
two possible channels, symmetric or anti- symmetric, which 
lead to the destructive interference of this coupling / e ff oc 
(Aj-A^/D — (A^A^/D. It is precisely the role of the mag- 
netic dipole-dipole interaction to split these channels, sup- 
pressing the perfect destructive interference, and leading to 

tm = w 2 - +t' 2 + ) -1 w) 2 > J % = ~b^ J ^ < 3) 

This Hamiltonian describes the flip-flop interaction between 
the 14 N nuclei leading to an exchange of the spin excitations. 
In Fig. 2(a), we present a scheme for the electron-mediated 
gate between two NV nuclei based on Eq. (3), referred as the 
nuclear XX gate. The initialization yields the state |l//b) = 
|0e) ® |<p n ) = |0,0) e ® |0, l) n , where electrons belong to the 
ground-state manifold of Fig. 1(b), and the dynamics of the 
spin excitation is determined by virtual electron spin-flip pro- 
cesses. In Fig. 2(c), we study numerically the accuracy of the 
effective Hamiltonian (3), which is compared to the exact evo- 
lution under the total Hamiltonian (l)-(2). One observes that 



the electron state remains in the ground-state, whereas there is 
a periodic exchange of the spin excitation between the nuclei. 
The remarkable agreement of both predictions justifies the va- 
lidity of the effective nuclear spin-spin Hamiltonian in Eq. (3). 
Unfortunately, the parameters in Table I yield a vanishingly- 
small coupling J^ « 0.1Hz, which is far too slow to produce 
any observable coherent coupling between the nuclei. Even if 
not of practical use, the above derivation gives a neat account 
of the mechanism of electron-mediated interactions, and will 
help us in understanding how to raise the interaction strength. 

A possibility to overcome this problem is to apply a mag- 
netic field, such that the Zeeman shift reduces D — » D — 
g e /lB#, thus enhancing J*^. Yet, one faces two important 
problems: i) In general, the axes of the NV centers are not 
aligned, and each electronic spin experiences a different Zee- 
man shift. For the large fields required, this inhomogeneity 
might exceed the dipolar coupling, and thus spoil the scheme. 
ii) The dephasing exerted by the environment would have a 
contribution that ruins the coherence of the interaction. We 
demonstrate below that there is a different approach that over- 
comes both problems simultaneously, and yet enhances the 
nuclear spin interaction: continuous microwave driving [18]. 

Effective driven interactions.- We discuss now the effects 
of a continuous microwave field that drives both the electronic 
and nuclear spins. The effect of the driving is two-fold: i) By 
addressing each NV center with different microwave fields, 
one can independently tune their frequencies so that they be- 
come resonant with a particular transition. This allows us to 
overcome the problems associated with both the inhomoge- 
neous broadening, and the different Zeeman shifts. More- 
over, this can be used for single addressing of NV's, especially 
when combined with magnetic gradients, ii) By tuning the mi- 
crowave frequency on resonance with the transition, one intro- 
duces a new energy scale that governs the system, namely the 
Rabi frequency. This parameter can be tuned by controlling 
the microwave power, allowing us to enhance / e ff . 

Let us consider the Zeeman effect associated to B = 30 G 
in Fig. 1(c). By setting the microwave frequencies to G) ej = 
Dj — ge/iB^j, C0nj = Pj — gn^Bj, one resonantly drives the 
transitions between the electronic and nuclear levels Mj = O 
— 1, Mj = O — 1. These driving terms can be written as 
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where the Rabi frequencies of the electronic and nuclear tran- 
sitions are £2 e >^n, and the electronic and nuclear Pauli matri- 
ces Gj, Tj. In the interaction picture with respect to #0,1 = 

ZjDj(S)) 2 -Pj(I Z j) 2 +g e LtBBjS Z j-gnVNBjIl one can neglect 
the rapidly oscillating terms associated to the transverse part 
of Zeeman shifts, and the hyperfine coupling. This rotating 
wave approximation is justified for the parameters shown in 
Table I. Additionally, we consider two NV centers with dif- 
ferent axes, which allows us to neglect the transverse part of 
the magnetic dipole coupling. For weak-enough driving, we 
arrive at the the total driven Hamiltonian 









Table I. Specific values of the coupling strengths 
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Figure 2. Nuclear quantum gate between two NV centers: Scheme for the initialization, evolution, and read-out of the effective J^ f 
interaction in (a), and J*f f interaction in (b). Here, Xj^ = exp(i0T^/2),7 / ^ = exp(i0Ty/2) represent finite pulses, and ^(^#) stand for the 
electron (nuclear) spin polarization, (c) Comparison between the effective dynamics under the Hamiltonian (3) and the exact time evolution 
under Eqs. (l)-(2) for the J*^ nuclear interaction. The expectation values represented correspond to the nuclear spin (I z .), and electronic spin 



(S z lo u 



(S\ + «S|). (d) Comparison between the exact (5) and effective (6) dynamics for the J*f f nuclear interaction, together with an echo 



scheme that allow us to get rid of the fast single-nuclei dynamics. We represent the nuclear expectation values 5(t? ) = (rf)t f — (tf)t - The 
dotted lines correspond to Ju = 0, where there is no interaction induced on the nuclei, (e) Performance of the ZZ-gate in the presence of 



different strengths of the electron dephasing noise bj = {5, 15,25,35,50,55}kHz, where the nuclear noise is Bj 
Ramsey decoherence times are roughly T 2e » {0.2, 0.07, 0.04, 0.03, 0.02, 0.018}ms, T 2n = 0.1 r 2e . 
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We stress that these approximations are justified by the pa- 
rameters in Table I, and supported by numerical simulations. 
We derive now the electron-mediated nuclear spin interac- 
tions starting from Eq. (5). We note that there is again a hier- 
archy in the couplings £l e ^> A 1 '- ^> Jn ^> £^ n , which leads to 
the clustering of energy levels shown in Fig. 1(d). By consid- 
ering the electron ground-state, the nuclear spins can interact 
through virtual electron spin-flips to the excited manifolds. In 
this driven regime, it is the longitudinal hyperfine coupling A 1 '- 
which induces such virtual transitions. A Schrieffer- Wolff- 
type calculation yields the nuclear Hamiltonian 



rrZZ_ tZZ Z T Z 



h£ftnTT 






eff " 



80, 



(s+4 



(6) 



where we considered the inhomogeneous broadening of the 
hyperfine couplings t, = 2[(A^) 2 - (A'|) 2 ]/£Vi 2 . This Hamil- 
tonian is an Ising magnetic interaction between the nuclear 
spins, which are additionally subjected to a transverse field 
due to the driving, and a longitudinal field due to the hyper- 
fine coupling. As advanced previously, we have been able 



to enhance the electron-mediated nuclear interaction, which 
becomes J*f f « 0. 1kHz for the parameters in Table I. Remark- 
ably, the strength of the nuclear spin interaction has increased 
by three orders of magnitude J*f f « lO 3 /^. 

In Fig. 2(b), we schematically describe the necessary in- 
gredients for the nuclear ZZ gate. The initialization consists 
of the electron (nuclear) spin polarization &*(^K), together 
with single- spin gates. & is obtained by the optical pump- 
ing cycle available for NV centers [5,6], whereas M is based 
on the techniques developed for the nuclear single-shot mea- 
surement [10], followed by the electron state-dependent fluo- 
rescence [5, 6]. Once polarized, |0,0) e |0,0) n , one applies 
unitary gates based on microwave pulses of different duration, 
Yj^ % = (I + it?) , Yj_z = (I — iZj) (also for the electron spin), 

which lead to \Yo) = | )e ® I — l-)n- The evolution of this 

state is dictated by the interaction-picture Hamiltonian (6), 
which leads to U t z * h = e-^i^e-^fe-^e+^o^i. Due to 
the longitudinal field, and the additional contributions of //o,i, 
the simple periodic exchange of the nuclear spin excitation 
shall be accompanied by fast oscillations. In order to observe 



neatly the effect of the interaction, one may perform a spin- 
echo sequence, such that the nuclear spins are inverted at half 
the gate time by a microwave pulse Xj y7t = ir? . In this case, 
the fast single-nuclei oscillations refocus after the spin-echo 
period tf, and one observes solely the effect of the interaction. 
In Fig. 2(d), we compare the effective description (6) to the 
Hamiltonian (5), which display a clear agreement. In partic- 
ular, when the echo period matches twice the ZZ-gate time 
k — 2^zz = 7t/2J*fi « 9ms, one finds a perfect excitation ex- 
change (Tf) : -1 -)► +1,(15) : + 1 "^ _1 - Note that for /12 = 0, 
this effect is completely absent. Finally, considering tf = t zz , 
and setting the echo pulse along the y-axis, the dynamics gen- 
erates a entangled nuclear state \\jfo) — I — )e ® | — y+y)n — ^ 
\Wf) = I — )e^(|-y+y)n + |+y- y )n)/v / 2. Once the gate has 
been performed, the nuclear operators (/j), (tJ ) must be mea- 
sured. Since the state-dependent fluorescence is particular to 
the electron spins, one should map the nuclear information 
onto the electrons, and then measure. This can be achieved in 
a quantum non-demolition fashion by using a microwave on a 
electron-spin transition conditioned to the nuclei [10]. 

Decoupling from decoherence.- So far, our discussion has 
focused on the idealized situation of isolated NV centers. 
However, every quantum system is inevitably coupled to an 
environment that degrades its coherence. This phenomenon, 
known as decoherence, must be seriously accounted for in 
solid-state materials, where the system-environment coupling 
is usually strong. In the particular case of NV centers, the 
major source of decoherence is the coupling to other impu- 
rity spins, such as single substitutional nitrogen electron spin 
(PI center) in type lb diamond [21], or 13 C isotopes in type 
Ha [8]. The microscopic description of the spin bath is an in- 
tricate many-body problem, and is a current subject of intense 
research. Here, we use a phenomenological model of the bath 
that yields a fluctuating magnetic field shifting the resonance 
frequencies. Due to the spin interactions, this effective field is 
modeled as a stochastic Ornstein-Uhlenbeck process [19, 22] 



#noise = Z{bj(t)S)+Bj(t)Pj 



(7) 



where bj(t),Bj(t) are random processes with autocorrelation 
(bj(t)bj(O)) = b 2 jQ- r J\ (Bj(t)Bj(O)} =B 2 j Q- R i t , where bj,BJ 
represent variance of the zero-mean gaussian distributions, 
and rj,Rj the inverse of their correlation times. In particular, 
the decoherence time of an electronic (nuclear) Ramsey ex- 
periment is given by 7^ = l/bj (72 n = 1/Bj). By considering 
the particular time-dependence of these stochastic processes, 
we numerically integrate the noisy dynamics, and average for 



N = 10 3 realizations of the random process. This allows us to 
study the effects of decoherence on the gate. 

For the slow XX gate (Fig. 2(a)), the limiting factor is 
the nuclear dephasing time, which can attain values of 72 n « 
10ms. Even for the purest samples, the coherence of the gate 
is completely lost much before the target time t xx « 4.5s is 
reached. Therefore, the performance of this gate is extremely 
poor. For the fast ZZ gate (Fig. 2(b)), not only the nuclear- 
spin dephasing, but also the electron-spin dephasing limit the 
gate accuracy. In the dressed-state basis (see Fig. 1(d)), the 
electron dephasing tries to induce a transition between the 
different manifolds, introducing additional noise in the nu- 
clei. However, due to the strong driving £l Q , these processes 
are partially suppressed. Additionally, a sufficiently strong 
nuclear driving, £l n ^> Bj, (bjA n JD^) 9 provides an additional 
decoupling mechanism that enhances further the gate perfor- 
mance. In Fig. 2(e), one observes the announced decoupling, 
since the gate performance at the target time t zz « 4.5ms is 
extremely good even for shorter electronic coherence times 
ranging from 7^ « 0.1ms to 7^ ~ 50/is. Due to the decou- 
pling mechanism, the gate accuracy will actually be limited by 
the decay times 7i e . Moreover, at this time scale, energy will 
be pumped into the system by the continuous driving. How- 
ever, note that this limitation can be overcome since T\ Q can 
be increased by orders of magnitude by cooling. Accordingly, 
one can achieve high fidelities. 

Let us finally note that the effective decoupling mechanism 
presented here can also be used to improve the electron-spin 
gates based on the direct dipole interaction [17]. In that case, 
the role of the microwave driving is to prolong dephasing 
times and to bring the two dressed electronic transitions to 
resonance to overcome the inhomogeneous broadening. 

Conclusions and outlook. - We have demonstrated the fea- 
sibility for engineering electron-mediated spin- spin interac- 
tions between the nuclei of two NV-centers. By continuous 
microwave driving, this scheme allow us to decouple from the 
electronic and nuclear dephasing sources, and increase the ef- 
fective interactions by three orders of magnitude magnitude 
thus achieving J^ f « 0.1 kHz for distances of existing pairs of 
NV-centers [17]. This scheme opens the possibility for the re- 
alization of quantum information processors, quantum simu- 
lators and quantum- sensors [23] on the basis of NV-centers in 
diamond. Finally, we would like to stress the generality of this 
scheme, which can be applied to other solid-state technologies 
that are candidates for quantum-information processing. 
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I. SUPPLEMENTARY MATERIAL 

In the following sections, we provide a detailed discussion of some technical aspects used to derive the main results presented 
above. In the first section, we derive the effective nuclear Hamiltonians in Eqs. (3)-(6) by means of the so-called Schrieffer- 
Wolff transformation. Finally, in the second section, we describe the phenomenological model used to study the effects of the 
decoherence on the nuclear-gate performance. 



1. Schrieffer- Wolff transformation and effective nuclear Hamiltonians 

In this section, we review the theory of quasi-degenerate perturbation theory as a tool to perform the so-called adiabatic 
elimination of the fast degrees of freedom in a quantum-mechanical system [1]. In particular, this technique allows us to derive 
the effective spin- spin Hamiltonians for the nuclei. 

Schrieffer-Wolff transformation. -We shall assume that the Hamiltonian is of the form H = Ho + XV, where XV is a weak 
perturbation to the unperturbed Hamiltonian //o = Ha+Hb. Here, Ha and Hb describe the slow and fast degrees of freedom 

H A \ja)=E a \afl), H B \ap) =E p \ap), (8) 

where \E a — E a /\ <C \Ep — Ep/\. Accordingly, the frequencies associated to the j3 — )■ f3 f transitions, (Oppi = Ep — Epi, are 
much larger than those of a — » a' , (O aa f = E a —E a >, and they can be adiabatically eliminated. By performing a canonical 
transformation U = e s , S^ = —S, one constructs an effective Hamiltonian // e ff = U^HU that only involves the slow degrees of 
freedom, PpH^Ppi = HJ? & 8ppt, where Pp = £ a \ a P)( a P\' The canonical transformation S, and the effective Hamiltonian H & ff, 
can be constructed to any order of the perturbative parameter A . This is usually known as the Schrieffer-Wolff transformation in 
condensed matter [2], although it also arises in a broader context [3,4]. To second order [1], one finds the following expression 

{a\H^\a') = {E a +Ep)8 aa , + <a/3|AV|a'/3> + \ £ (afS\lV\a"fS")(a"f5"\XV\a f {S)( j— + - -j- Y (9) 

Since we have fast and slow components copp// ^> (O aa n, we may expand {(O aa n + C0pp//)~ l « 0)7l„(l — CQ aa "/cQpp//) in Eq. (9), 
which allows us to get an effective Hamiltonian that only acts on the slow degrees of freedom 

pit ^jSjS" 

In this manuscript, we make use of this canonical transformation to adiabatically eliminate the fast degrees of freedom of the 
electrons localized around a Nitrogen- vacancy impurity in diamond, and obtain an effective Hamiltonian for the slow 14 N nuclei. 
Effective nuclear Hamiltonian for the static regime.- According to the above discussion, one should first identify the fast and 
slow degrees of freedom of the Hamiltonian in Eqs. (l)-(2). Regarding the parameters in Table I, together with the energy levels 
in Fig. 1(b), it is clear that the fast degrees of freedom correspond to electronic spin flips j3 = {m\ , m^\, whereas the slow degrees 
of freedom involve the nuclei a = {Mi,M2}. Besides, the weak perturbation corresponds to the hyperfine interaction, which 
couples different electronic manifolds. We are interested in the dynamics within the electronic ground-state manifold j3 = {0, 0}, 
where the Schrieffer-Wolff transformation allow us to get the following Hamiltonian 

(a|// e ° f ?|a') = -E^(/)) 2 + ^ E {00,aXHf^\p''a'W'cc''XHt n) \00,a')( \ + * Y (11) 



where the energy zero is shifted to Eq = £y |(P y - — Dj), and one has to sum over all the electronic states H/3")} on the first- 
excited manifold in Fig. 1(b). There are thus two types of channels of virtual electron spin flips, and these correspond to 
the symmetric |S±) = (|0,±l) e + | ± l,0) e )/v / 2, and anti- symmetric |A±) = (|0,±l) e - | ± l,0) e )/v / 2 combinations. Due to 
the dipole-dipole interaction, these channels are split in energies by E s —E a = Z/12, which is the key ingredient that allow us 
to obtain a non- vanishing interaction between the nuclear spins. By expanding to leading order for fi" — S±,A± 9 we obtain 
(&W' + co oof5")~ l ~ ^oob" ~ D~ l (l ^fJn/D) for the symmetric/anti- symmetric states. By resuming the expression above, one 
finally arrives to the effective nuclear Hamiltonian, which reads 

ffeff = -L^) 2 +Wl + '2~+A~' 2 + ), (12) 

J 

where J*f f = lA^A^Jn/D 2 , and we have a negligible local energy shift for (Aj~) 2 /D <C D. This is precisely the spin-spin 
Hamiltonian H^ described in Eq. (3), which leads to the exchange of spin excitations between the nuclei. Let us finally remark 
that in this derivation we have assumed that there is no inhomogeneous broadening D\ = D2, which may modify the results. 
However, since the effective interactions are so weak J^ « 0.1 Hz, there is no point of being more rigorous at this point. In the 
next section, we shall treat the possible effects of inhomogeneous broadening for the driven interactions in detail. 

Effective nuclear Hamiltonian for the driven regime.- In this part of the Appendix, we discuss the RWA leading to Eq. (5), 
and the Schrieffer- Wolff transformation to the effective nuclear Hamiltonian in Eq. (6). 

a) Driven Hamiltonian.- Let us rewrite the total driven Hamiltonian in Eqs. (1),(2), and (4) as follows H = #0,1 +#0,2, where 

#0,1 ^^(S^O^^eMB^COSe^. -Pj(fy 2 -g n llsBCQs0jFj, 

3 < x , ^ (13) 

#0,2 = £ ifismfyteeMBe-^St - g^^J if + H.c.) +E 7 -#f n) +H[ Q 2 e) +H d (t), 

j 

where we again shifted the energy zero to Eq = £y | (Pj —Dj), and we have introduced the relative orientation of the NV centers 
(0/, <py) with respect to the applied magnetic fields. In the interaction picture #0,2 (0 = e l//o ^ifo,2e _l//oi/ , one can neglect rapidly 
oscillating terms by a rotating wave approximation (RWA), which leads us to 



2 0,2 



* £ (> e CTj + i<2 n Tj) + 2/i2^S| + £A^, (14) 



where we define Bj =Bcos Bj. This RWA is justified when Dj ^> g^ji^Bj ^> Q e ,Aj-, P ; ^> gnjl^B ^> £2 n , which is clearly fulfilled 
for the parameters shown in Table I. Additionally, the inhomogeneous splitting allows us to neglect the transverse part of the 
magnetic dipole coupling between the electron spins J\2 <C g Q jlQ\B\ — #2!- According to Table I, one can neglect these terms 
Aj/Dj ~ 1(T 3 , gzUBB/Dj - 1(T 2 , a e /g e fi B B ~ 10" 1 , g n ^B/Pj ~ 10" 3 , £l n /g n il N B - 10" 1 . Besides, for two NV centers 
oriented along different axes, Q\ — 62 ~ &{k), we can neglect the transverse dipole-dipole coupling /12/^eMB^i — Bl\ ~ 10 -3 . 
In order to confirm the validity of these approximations, we must compare the dynamics of both Hamiltonians, //o,2(0>#3T- 
Let us note that the full time-dependent Hamiltonian #0,2 (0 = t m ^ lt HQ^~ m ^ lt contains very different time-scales, ranging 
from ns to ms. To reproduce the dynamics faithfully, one sets the numerical integration time-step to the smallest time-scale, ns. 
For such a small time-step, prohibitively large integration times are required in order to reach the ms-regime where the nuclear 
spin-spin interaction effects become visible. Nonetheless, to test the accuracy of the RWA, it suffices to study t G [0,2k/ A j], 

which lies in the /is-range. In Fig. 3(a), we compare both predictions numerically for the initial state \\jfo) — I )e ® I — l-)n, 

namely 

<t?(0U = (Vole^^e-^TVo), (^(r)) exact = Me^^^e"^^^. (15) 

Since there is no refocusing echo pulse, the nuclear spin dynamics should be dominated by the Rabi oscillations caused by 
the term — \A l -T z - in Eq. (6). As observed in Fig. 3(a), these neat jis Rabi flops display a perfect agreement between the exact 
Hamiltonian and the RWA approximation. 

b) Effective nuclear Hamiltonian: According to the preceding discussion, we shall consider directly the RWA Hamiltonian in 
Eq. (5), where the set of Pauli matrices is defined as follows 

^=|0;>e(0;|e-|-l;>e(-l;|e, ^=|-l y -)e(0 7 -| e +H.C. 

T z j=\oMOj\ a -\-iM-ij\ a ,T!=\-iMOj\ n +n.c. 

and the sub-indexes | ) e , | ) n indicate the electronic or nuclear origin of the spin state. In this two-level approximation, the 
spin-1 operators become S z - = | ly) e (ly| e + \ (of — ^2)* and I z = \lj) n (lj\ n + ^(t z - —I2). Accordingly, the electronic and nuclear 
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Figure 3. Rotating wave approximation and free induction decay, (a) Accuracy of the rotating wave approximation: Dynamics of the 
expectation values {z x (O)rwa (squares, circles), as solved by numerical exponentiation, and {z x (O)exact (solid lines), as solved numerically by 
a fourth order Runge-Kutta method, (b) Free induction decay due to magnetic noise: Numerical simulation of the stochastic dynamics of the 
free induction decay ("C x (t)} (circles), for the initial state I^Pq) = (|0) e + | — l) e )/V2, after averaging for N& = 5 • 10 3 trajectories of the random 
process (25). The red line corresponds to a gaussian fit (^ x (t))m ^ ex P( — i^lt^ 2 )' where b^_ = 1.09kHz. 



levels rrij = 1,M/ = 1 decouple, and one may write the following driven pseudo spin- 1/2 Hamiltonian 

if =L(^e^ + ^ n T-)+^ 12 (CTf-I 2 )(CT|-I 2 ), Hi=E^A^-I 2 )(T?-l2), 



(17) 



Once the pseudospin-1/2 Hamiltonian in Eq. (17) has been derived, we adiabatically eliminate the fast electronic degrees of 
freedom from the slow nuclear dynamics by a Schrieffer- Wolff transformation. We shall make use of Eq. (10), where we identify 
j3 = { — ,—} as the lowest-energy manifold (see Fig. 1(d)), j3 7/ = S, A as the symmetric/anti- symmetric excited manifolds 
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(18) 



where we have introduced |±) e = (|0) e ± | — l) e )/\/2, and a parameter £ <C 1 quantifying the inhomogeneous broadening A'J 7^ 
Aj. These levels are split in energies by 2J(t;), where 



To second order in the hyperfine coupling, we are able to derive 



(19) 



fl*=B«-*y + ^>*J) - i6^T< — II>W*/l s XSII>!t^^l — > - i^-< — II>WT/l A XA|I>!tc^^| — >, (20) 
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where we have introduced the following energy differences 
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(21) 



By computing the corresponding matrix elements, together with a Taylor expansion for £l e ^> A'j, /12, and 1 ^> §, we find the 
following expression for the followig nuclear spin Hamiltonian, which is precisely Eq. (6) in the main text, 



n; e = j^\x\ +lfi n T* - l -A>fi, j% =-§^( J -£ -r 1% !• 



oi2 e Y"e 



(22) 



2. Decoherence and effective decoupling by continuous microwave driving 



In order to perform quantum-information tasks in a solid-state device, the effects of the system-environment coupling must be 
carefully addressed. In contrast to cold-atom platforms, the environment in a solid is rather complex since the spins may couple 
to a wide variety of excitations. In the case of NV centers, whose energy levels lie deep in the band gap of diamond, the major 



source of noise is the coupling to the spins of different impurities, rather than to electronic or vibronic excitations. Accordingly, 
one should consider the effects of a spin bath on the coherent features of the electron/nuclear spin of the NV center. 

Phenomenological magnetic noise model - The problem of a central spin coupled to an ensemble of bath spins has been 
studied since the early days of nuclear magnetic resonance [5], and depending on the particular nature of the spin bath can be an 
intricate many-body problem. For type lb diamond, the bath consists of the electronic spins of 14 N impurities, the so-called PI 
centers, randomly distributed through the sample. The dipolar coupling of the PI centers to the NV electron spin gives rise to 
a pure dephasing which can be treated by mean-field theories [21]. Conversely, for ultrapure type Ha diamond, it is the nuclear 
spin of 13 C isotopes which yields the dephasing of the NV center via the hyperfine electron-nuclei coupling [8]. Interestingly, 
the correlations of this nuclear-spin environment must be accounted in order to reproduce the short-time dynamics of the system. 
In this work, we follow a phenomenological approach rather than a microscopic one, where the magnetic noise is modeled 
by a random fluctuation of the resonance frequencies associated to the electron/nuclear spins. This model captures the whole 
dynamics of type lb diamond [19], and is expected to describe faithfully the long-time dynamics of type Ha diamond (tf ~ms), 
where the non-Markovian aspects of the environment should not have an important effect. 

The flip-flop interactions between the NV and the bath spins can be safely neglected due to their utterly different energy 
scales. Therefore, we consider that the collective effect of the spin bath is to shift the resonance energies of the NV center (i.e. 
pure dephasing), which can be modeled by an effective local magnetic field. In order to account for the spin-bath dynamics, this 
magnetic field is treated as a stochastic process [19]. Hence, the NV centers are described by a stochastic Hamiltonian 

H({bj(t),Bj(t)}) = £ (#f +Hf +tfj e " n) ) +H^ e +H d (t)+H noise , H noise = £ (b j (t)S c j +B j {t)lf) , (23) 

j J 

where bj(t),Bj(t) are the random processes describing the fluctuating resonance frequencies. Whereas the electron and nuclear 
spins of the NV center follow the unitary dynamics dictated by the Schrodinger equation i<9 / | v F) = if ({&,■(*), 5/(7)}) |*P), the 
fluctuating fields evolve according to the so-called Langevin equation d t X(t) =A(X(t),t)X(t) + ^/D(X(t)^t)T(t), where X(t) = 
{bj(t) : Bj(t)} is the random process, A(x,t) 9 and D(x,t) > 0, are smooth functions, and T(t) is a Gaussian zero-mean noise (i.e. 
(r(f )) = 0, (T(t)r(t')) = S(t — t')) [6]. Due to the large number of spins conforming the bath, these random processes can be 
argued to be Gaussian by means of the central limit theorem. Besides, when the back-action of the system is small [19], or we 
are interested in the long-time dynamics, the process can be treated as Markovian and stationary. In this case, one obtains a 
Langevin equation with A(x,t) = —x/t, and D(x,t) = c, where % is the relaxation time of the process, and c the diffusion speed. 
The solution to this equation can be obtained explicitly, and is known as the Ornstein-Uhlenbeck (OU) process [7]. In the case 
of zero-mean random magnetic fields, the auto-covariances are 

(bj(t)bj(O)) = by J', (Bj(t)Bj(O)) = By R J', (24) 

where b] = Ic^.T^.^B 2 - = ic».TH. represent the variances of the zero-mean gaussian distributions, and r 7 - = \/%.,Rj = 1/t». 
the inverse of their relaxation times. These auto-correlations lead to a Lorentzian spectral density, which contains a white-noise 
region at low frequencies, and l// 2 -noise region at larger frequencies. Interestingly enough, the dynamics of the OU process 
can be given explicitly [7] , and numerical integration of the Langeving equation is not required. In fact, for any discretization 
dt > 0, one finds the following exact update formula 



X^ + d0=X(0e-^ T + Wy(l-e- 2d ^K (25) 

where n is a zero-mean unit- variance gaussian random variable which is time uncorrelated. In order to solve the whole stochastic 
quantum dynamics in Eq. (24), we discretize the time interval in M time-steps, t m = mdt G [0,ff], where dt = tf/M, and obtain 
the different values of the fluctuating magnetic fields bj(t m ),Bj(t m ) by employing the above formula (25). Then, we integrate 
numerically the stochastic Hamiltonian for the particular sampling of the random process s = {bj(t m ),Bj(t m )}, and recover the 
expectation values (tj) s . By repeating this procedure for N{ t ^> 1, one can perform the statistical average over the stochastic 

noise, (l x - ) = ^- £ 5 (?y ) s * an d thus study the effects of the decoherence. 

To illustrate the physics of this phenomenological model, let us consider the simpler situation of a single NV electron spin. 
We consider the decoherence of a Ramsey experiment, where the initial state corresponds to |\Fo) = (|0) e + | — l)e)/\/2, and 
we measure the so-called free induction decay (FID) due to the noise after a certain time t, (o x (t)). In Fig. 3(b), we represent 
the time evolution of the FID derived from the numerical solution of the stochastic Hamiltonian H(b(t)) =b(t)S z , where b{f) is 
a OU process with b = 1kHz, and we have averaged over N{ t = 5 • 10 3 samplings of the random process. Due to the magnetic 
noise, the free induction decay follows a gaussian decay law (T x (t)) <* exp( — ^b 2 t 2 ), which allows us to identify the dephasing 
time as 7^ = 1/b & 1ms. Therefore, we observe how the phenomenological noise model allows us to study the decoherence 
effects for different dephasing rates, which has been used in the main text (Fig. 2(e)). 

Effective decoupling mechanisms.- An advantage of the phenomenological noise models is that they allow a neat under- 
standing of the effects of decoherence, together with possible strategies to overcome them. In the particular case of the driven 



Hamiltonian (5), it is easy to observe that the fluctuation of the electronic resonance frequencies, H = Y,jbj(t)Sp tries to induce 
transitions between the energy manifolds of Fig. 1(d), namely |+ 7 -)e ** |— j)e- However, since these states now have a huge en- 
ergy difference given by the Rabi frequency of the driving £l e , these transitions are non-resonant and thus partially suppressed. 
In fact, the electron magnetic noise can only couple to the nuclei via second order processes. The leading order contribution 
comes from the coupling to the hyperfine channel, and gives rise to H e ff « £ / -(£ / -A l '-/Q e )Ty, which is partially suppressed for 
the regime considered in this work bj,Aj <C £2 e - Now, one has to compare this new term to the nuclear driving, and since 
(bjA u j/Q, G ) < (l n , we get an additional decoupling mechanism. Qualitatively, one can argue that the effects of the noise give 

rise to a small second-order fluctuation of the nuclear driving H e ff « E 7 -n n (l + ^(bjA n j/D^) 2 /Q^)tj. In this expression, one 
observes the two-fold role of the microwave driving £l Q . On the one hand, £l Q must be small enough so as to increase the effective 
nuclear interaction. On the other hand, £2 e must be big enough so as to provide an effective decoupling from the electronic noise. 
Therefore, one must find a compromise between the two, such as that presented for the parameters in Table I. With respect to 
the additional decoupling due to the driving of the nuclei, £l n must be as big as possible. By increasing the external magnetic 
fields beyond B « 500G, where the levels rrij = 0,-1 become degenerate, one could raise the nuclear driving strength, and thus 
increment the efficiency of the decoupling. 
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